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DATA PROCESSING, PARTICULARLY IN COMMUNICATIONS SYS- 
TEMS 

Field of the Invention 

This invention relates to methods for solving systems of linear equations and, in 
particular to increasing speed with which such solutions can be obtained by using a 
computer system. This invention also relates to a computer system for solving systems 
of linear equations, the computer system without multipliers and dividers. This inven- 
tion also relates to applications of the method and computer system in communications 
and signal processing, specifically to the applications which require solving systems, of 
linear equations arising from the linear least squares (LS) problem. More specifically, 
this invention also relates to multiuser receivers and adaptive filtering. 
Background of the Invention 

Linear equations occur frequently in all branches of sciences and engineering, and 
effective methods are needed for solving them. Furthermore, many science and engi- 
neering problems lead not to a single equation, but to a system of equations. Classical 
methods of solving these systems can be divided into two categories: direct methods 
and iterative methods. 

Direct methods attempt to produce an exact solution by using a finite number of 
operations. A problem with direct methods, however, is that the number of operations 
required is large, which makes the methods slow and, also, sensitive to truncation er- 
rors. Furthermore, direct methods often fail on ill-conditioned matrices. Also, due to 
truncation errors and complex operations required such as division and multiplication, 
direct methods are not well suited for implementation on fixed-point computer systems. 
Instead, direct methods are implemented preferably on floating-point computer systems. 
However, floating-point computer systems are slower, more complicated, and more ex- 
pensive for implementation compared to fixed-point computer systems. 

Iterative methods solve a system of equations by repeated refinements of an initial 
approximation until the result is acceptably close to the solution. Each of the iterations 
is based on the result of the previous one, and in theory, is supposed to improve the pre- 
vious approximation. Generally, iterative methods produce an approximate solution of 
desired accuracy by yielding a sequence of solutions, which converges to the exact solu- 
tion as the number of iterations tends to infinity. 
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In solving systems of equations, especially when the number of equations is large, 
it is desirable to use a computer system. The word length of the system has a direct 
bearing on accuracy, and the likelihood of truncation errors increases with the number 
of operations required for a solution. For this reason, iterative methods are often pre- 
5 ferred over direct methods because a solution can be arrived at with fewer operations. 
Yet, existing iterative methods do not adequately minimise the number of operations 
required to reach a solution. Another disadvantage of existing iterative methods is that 
they are not able to avoid truncation errors even if the number of operations is small. 

When solving a system of linear equations with a computer system, another con- 
10 sideration is hardware efficiency. For example, it is often desirable to implement itera- 
tive methods in hardware, such as FPGA (Field-Programmable Gate Array) or ASIC 
(Application-Specific Integrated Circuit) platforms. One way to improve efficiency is to 
exclude operations, which are complex for hardware implementation. One skilled in the 
field knows that multipliers and dividers occupy areas on FPGA or ASIC, which are one 
15 order greater than that of adders. Existing iterative methods exploit multiplication 
and/or division operations. Therefore, existing iterative methods result in non-efficient 
hardware implementation. 

Another way to improve efficiency is the use of a computer system with fixed- 
point operations. However, existing iterative methods, when implemented in such sys- 

20" terns, lead to low accuracy of the solution: Furthermore, truncation errors can lead to 

divergence of the existing iterative methods. 

Another shortcoming of existing iterative methods is that the number of iterations 
is not decreased when the solution is sparse; here, the sparse solution means that only a 
relatively small number of elements of the solution vector with respect to the vector 
25 length are not zero or have high amplitudes. Many sciences and engineering problems 
deal with sparse solutions. For example, in communications when an impulse response 
of a channel with multipath propagation is estimated, the number of valuable paths can 
be as small as between 1 and 6, while the number of resolvable delays, defining the size 
of a system of linear equations, can be as large as several hundreds. 
30 In communications and signal processing, there are many applications where solv- 

ing systems of linear equations is required. This includes linear equalisation, for exam- 
ple, zero-forcing and MMSE (minimum mean squares error) equalisation. This also in- 
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cludes measurement of an impulse response of a propagation channel, for example, a 
multipath channel, a room reverberation, and others. This also includes linear multiuser 
detection, for example, decorrelating and MMSE multiuser detection. This also includes 
speech coding, for example, linear predictive coding (LPC). In these and many other 
5 applications a system of linear equations may arise from the least squares (LS) problem. 

CDMA (Code Division Multiple Access) is a multiple access method based on a 
spread spectrum technique. CDMA has many advantages over earlier developed FDMA 
(Frequency Division Multiple Access) and TDMA (Time Division Multiple Access) 
methods, such as simplicity of frequency planning and spectrum efficiency. In a CDMA 

10 method, a narrow-band data signal of a user is multiplied to a relatively broad band by a 
spreading code having a much broader band than the data signal. All users transmit si- 
multaneously on the same frequency band. On each connection between a base station 
and a mobile user, a different spreading code is used, and the signals of the users can be 
distinguished from one another in the receivers on the basis of the spreading codes of 

15 the users. 

The object of the receiver is to detect the signal of the desired user among a plural 
number of interfering signals. In the receiver, the data signal is restored to the original 
band by multiplying the received signal, in correlators or matched filters, by the same 
spreading code as was used in the transmission step. Ideally, the signals of the other us- 

20 ers are not restored to the narrow band: In practice, the spreading codes may partly cor- 
relate, and in such instances the signals of the other users make it more difficult to de- 
tect the desired signal, by distorting the received signal. This interference caused by the 
users to one another is called multiuser interference. 

Multiuser receivers are effective means for cancelling multiuser interference and 

25 such cancelling is known to improve spectral efficiency of multi-user communications 
networks. 

One receiver solution operating in an environment with multiuser interference is a 
so-called decorrelating receiver, which cancels multiuser interference by multiplying 
matched-filter outputs by the inverse matrix of the cross-correlation matrix of the 
30 spreading codes. A decorrelating receiver is discussed in greater details in R. Lupas, S. 
Verdu: "Linear multiuser detectors for synchronous code-division multiple access chan- 
nels", IEEE Transactions on Information Theory, vol. IT-35, p.123^136, January 1989. 
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In a channel without multipath propagation the cross-correlation matrix can be precom- 
piled and inverted, and the inverse matrix can be stored in the decorrelatirig receiver; 
this enables the receiver to avoid calculating the inverse matrix and, thus, simplify the 
receiver. However, in known decorrelating receivers, multiplying by the inverse matrix 
5 is complicated because it requires multiplication operations. Another disadvantage of 
the decorrelating receiver is that in channels with multipath propagation, i.e., in more 
realistic channels, the cross-correlation matrix depends on multipath characteristics, 
therefore the inverse matrix cannot be precomputed and matrix inverse operation is re- 
quired which is very complicated for real-time implementation, especially for a large 

10 number of users. Another disadvantage of known decorrelating receivers is the effect of 
noise enhancement leading to increasing detection errors. 

One skilled in the field knows that a better detection performance is achieved in a 
so-called MMSE (Minimum Mean Squares Error) multiuser receiver. A disadvantage of 
the MMSE receiver is that it uses a modified cross-correlation matrix which depends on 

15 the noise power and whose inversion cannot be precomputed. Another disadvantage is 
that, in known MMSE receivers, multiplying by the inverse matrix is complicated be- 
cause it requires multiplication operations. 

Many multiuser receivers (such, for example, as decorrelating receiver, MMSE re- 
ceiver, and interference cancellation receivers) are implemented on the basis of solving 

20 ^ linear systems of equations; this is known, for example, from the paper by Peng H. Tan 
and Lars K. Rasmussen, "Linear interference cancellation in CDMA based on iterative 
techniques for linear equation systems", published in the journal IEEE Transactions on 
Communications, vol. 48, No. 12, Dec. 2000, pages 2099-2108. In multiuser receivers, 
direct methods for solving systems of linear equations are complicated for implementa- 

25 tion. Existing iterative techniques, such as, for example, methods with the steepest de- 
scent iteration, the conjugate gradient iteration, the Jacobi iteration, the Gauss-Seidel 
iteration and others are yet complicated for implementation. Therefore, multiuser re- 
ceivers are still limited in use, which does not improve spectral efficiency of multi-user 
communications network, for example, the third generation mobile communication sys- 

30 tern UMTS, and the like. When used in equipment of a mobile user, for example, in 
UMTS downlink, multiuser receivers complicated for implementation result in a high 
power consumption, which shortens the battery life of the equipment. 
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Adaptive digital filtering is used in many areas of signal processing and communi- 
cations. Optimal adaptive filtering is based on solving systems of linear equations; for 
example, this is known from the paper by John F. Doherty, "Channel equalization as a 
regularized inverse problem" published in the book "The digital signal processing 
5 handbook", CRC Press, 1998, page 31-3. Existing direct methods for solving systems of 
linear equations are complicated for implementation. In practice, iterative methods such 
as recursive least squares (RLS) and least mean square (LMS) methods are most suited 
for real-time implementation. The RLS method possesses a fast convergence, however, 
it is computationally unstable and has a high complexity; in particular, it requires multi- 
10 plication and division operations. The LMS method can be implemented with no divi- 
sion, but it still requires multiplications. Another disadvantage of the LMS method is 
that it possesses a low convergence speed, which significantly limits its applications. 

Since existing iterative methods exploit multiplication operations they are compli- 
cated for hardware implementation, this limits a sampling frequency at which adaptive 
15 digital filters can operate and applications where such filters can be used. For example, 
it is a practical problem of implementing adaptive filters operating at sampling frequen- 
cies as high as 100 MHz and higher. Another disadvantage of existing adaptive filters 
based on iterative methods is that elements of the solution vector are truncated before 
using them as filter coefficients; this results in a truncation error of adaptive filtering. 
20" Another disadvantage of the RLS and LMS methods is that after the convergence there 
is a final error with respect to an accurate solution; the error is determined, in particular 
by a constant step-size parameter. 
Summary of the Invention 

One aspect of the invention is a method for solving systems of linear equations. 
25 Elements of a solution vector are represented as fixed-point binary words. The method 
comprises the steps of initialising a solution vector and an auxiliary vector and bit-wise 
iterations performing passes through elements of the solution vector, updating elements 
of the solution and auxiliary vectors in the passes, and repeating the passes until a re- 
quired accuracy of the solution is achieved. The method solves a system of linear equa- 
30 tions by starting the bit- wise iterations with the most significant bit of the words and 
proceeds to a next (less significant) bit when the solution achieved for the current bit is 
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not changed in the last pass. The solution is obtained with a computer system pro- 
grammed to perform a sequence of operations. 

Another aspect of the invention is a computer system for solving systems of linear 
equations, the computer system performing no multiplication or division operations. 
5 The computer system comprises: a host processor producing itself or receiving from 
other devices parameter signals representing elements of a coefficient matrix and right- 
side vector of the system of linear equations and transmitting to other devices parameter 

c 

signals representing elements of a solution vector; a host bus coupled to the host proces- 
sor; an internal bus; a first means for updating and storing elements of the solution vec- 

10 tor, the first means coupled to the host and internal buses; a second means for storing 
elements of the coefficient matrix, the second means coupled to the host and internal 
buses; a third means for determining successful iterations and preferable updates, the 
third means coupled to the internal bus and the second means; a fourth means for updat- 
ing and storing elements of an auxiliary vector, the fourth means coupled to the host 

15 bus, internal bus, and the second means. 

Another aspect of the invention is a multiuser receiving method in a data transmis- 
sion system in which code division multiple access involving multiuser interference 
among respective signals, each signal representing a succession of data signals trans- 
lated into bits and transmitted at a rate of a plurality of chips per bit, spread by a respec- 

'20 " tive spreading code, is applied for detecting a particular data signal, from among a plu- 
rality of data signals, said method comprising: filtering matched with the spreading 
codes and applied to the received signal to obtain respective output signals; transform- 
ing the matched-filter output signals by solving a system of linear equations of the kind 
Rh=p where R is a NxN cross-correlation matrix of the spreading codes, P is a Nxl vec- 

25 tor grouping the matched-filter output signals, h is the Nxl solution vector representing 
the transformed signals, and TV is a number of used spreading codes, and subjecting the 
transformed signals to obtain an estimate of the data signal; wherein solving of the sys- 
tem of linear equations comprises the steps of: representing elements of the solution 
vector as fixed-point binary words each consisting of at least one bit; initialising the so- 

30 lution vector and an auxiliary vector; performing, for each bit representing the binary 
words, bit-wise iterations comprising the steps of: performing passes through all ele- 
ments of the solution vector; updating elements of the solution vector in the passes; up- 
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dating elements of the auxiliary vector in the passes; repeating the passes until a finish- 
ing condition is fulfilled; and performing the bit-wise iterations until a stopping condi- 
tion is fulfilled. 

Another aspect of the invention is a multiuser receiver performing matched filter- 
ing of a received signal with all user-specific spreading codes and transforming result of 
the matched filtering to a solution vector whose values are represented by a solution of a 
system of linear equations; the multiuser receiver comprising means to solve the system 
of linear equations without performing multiplication or division operations comprising: 
means for matched filtering, a computer system for generating a solution vector h in an 
equation of the kind Rh=p, where R is the NxN cross-correlation matrix of spreading 
codes, P is the Nxl matched-filter vector, and h is the Nxl solution vector whose ele- 
ments are to be calculated, and means for estimating data signals from the output signals 
of the computer system. 

Another aspect of the invention provides an adaptive filter having a transfer func- 
tion determined by a plurality of coefficients whose values are represented by a solution 
of a system of linear equations; the adaptive filter comprising means to solve the system 
of linear equations without performing multiplication or division operations, compris- 
ing: means for calculating autocorrelation and cross-correlation coefficients, means for 
generating a solution vector h in an equation of the kind Rh=P, where R is the NxN co- 
efficient autocorrelation matrix, p is the Nxl cross-correlation vector, and h is the Nxl 
solution vector whose elements are to be calculated. 

The invention also provides digital signal processing apparatus comprising circuit 
means adapted to receive a signal and to perform an operation upon a said signal in such 
a way as to alter a significant parameter thereof, the operation being conducted in re- 
sponse to a method according to the invention for solving systems of linear equations. 

The invention further provides digital signal processing apparatus, comprising cir- 
cuit means adapted to receive a signal in the form of a plurality of signal paths, repre- 
sented by a plurality of simultaneous linear equations as a function of a significant pa- 
rameter of the said signal, calculating means to calculate the said coefficients and to 
provide an output in which the signals have the said significant parameter altered ac- 
cording to the calculated coefficients, the calculating means operating under the control 
of a method according to the invention for solving systems of linear equations. Typical 
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examples of digital signal processing apparatus are: an adaptive transversal filter; an 
echo canceller, a steerable antenna array; a computed tomography apparatus; apparatus 
for correcting multipath distortion, such as in RAKE receivers; a speech processor, such 
as a linear predictive coder; apparatus for modifying room acoustics; or digital image 
5 processing apparatus, such as a television image processor. It will be appreciated, how- 
ever, that the invention is not limited to such examples. 

A technical advantage of the invention is that a system of linear equations may be 
solved with a minimum of operations. A further advantage of the invention is that solu- 
tion may be obtained without multiplication and division operations, enabling , simple 
10 hardware implementation. A further advantage of the invention is that solution may be 
obtained with higher accuracy. Another advantage of the invention when applied to 
multiuser detection is that spectral efficiency of communication networks is increased 
and power consumption of multiuser receivers is reduced. Another advantage of the in- 
vention when applied to adaptive filtering is that the convergence speed increases and 
15 truncation errors are eliminated. 

In practical terms, the avoidance of multiplication steps in the algorithm translates 
into there being no necessity for multipliers in a circuit or chip designed to implement 
the algorithm. When it is appreciated that a multiplier demands a complexity that is an 
order of magnitude greater than that of an adder, it is immediately evident that substan- 
20 tial savings can be made by the use of the algorithm. This leads to consequential bene- 
fits in terms of speed, compactness, reduced power consumption, and improved port- 
ability of devices incorporating the invention. 
Brief Description of the Drawings 

The invention will now be described with reference to the drawings, which illus- 
25 trate exemplary embodiments of the invention, and in which: 

Figure 1 is a flow diagram implementing the method for real-valued data; 

Figure 2 is a flow diagram implementing the method for real-valued data and an 
auxiliary delta array; 

Figure 3 is a flow diagram implementing the method for real-valued data, an auxil- 
30 iary delta array, and a check for the solution amplitude range; 

Figure 4 is a flow diagram showing application of the method for complex-valued 

data; 
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Figure 5 is a flow chart explaining the steps in searching for the minimum; 
Figure 6 is a flow diagram showing application of the method for complex-valued 
data and an auxiliary delta array; 

Figure 7 is an example of a data processing system implementing the method; 

Figure 8 is a general scheme of the computer system according to the invention; 

Figure 9 represents an R-memory block; 

Figure 10 represents an h-updates block; 

Figure 1 1 represents a Q-updates block; 

Figure 12 represents a minimisation block; 

Figure 1 3 represents a multiuser receiver; 

Figure 14 represents a detector of a multiuser receiver; 

Figure 1 5 represents an adaptive filter; 

Figure 1 6 represents a block for calculating filter coefficients in an adaptive filter. 
Detailed Description of the Illustrated Embodiments 

The invention describes a method for solving systems of linear equations, a com- 
puter system for solving systems of linear equations without multiplication or division 
operations, a multiuser receiving method, a multiuser receiver, and an adaptive filter. 

Prior to description of the preferred embodiments, a description of the least squares 
problem and coordinate descent optimisation will be made so as to help better under- 
"standihg of "the present invention. 

The linear least squares (LS) problem deals with minimisation of the function 

J(h) =| Zh - d | 2 

with respect to an unknown Nxl vector h, where Z is an MxN matrix, and d is a 
Mxl vector. This problem is known to be equivalent to solving a system of linear equa- 
tions (normal equations): 

Rh = p 

where R is an NxN coefficient matrix 
R = Z T Z 

P is an Nxl right-side vector 
p = Z T d 

and Q T denotes matrix transpose. The function J(h) can be written as 
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J(h) = i £ £ R(m, n)h(m)h(n) - £ P(m)h{m) 

where /2(m; is 772 th element of the vector h, J3(m) is /77 th element of the vector p, and 
R(m,n) is {m.nf^ element of the matrix R. 

In practical circumstances, elements of the solution vector h are limited in ampli- 
tude. Therefore, minimising the function J(h) with respect to an A^-dimensional vector h 
is considered under constraints 

h e U = {[*(!),..., h(m),... t h(N)] :| h(m) \<H,m = 1,..., A^} 

where //>0 is a known number. 

A number of computerised techniques are known for minimising a function J(h); 
coordinate descent optimisation is one of them. The coordinate descent optimisation can 
be described as follows. Let = (0,...,1,...,0) be a unity basis vector, whose z ,th coordi- 
nate is 1 and the others are 0. Let h 0 be an initial value of the vector h and czo be a posi- 
tive number. Let us also assume that we know a vector h k and a number Gfc for some 
k>0. Denote 

where [k/N] denotes the integer part of the number k/N; this means that 

Pi = e i? P2 = e 2' * *' Pn = 

Pn+] = e i? Pn+2 = e 2 5 —5 P2A^ = C AM 

Let us calculate the function J(h) at the point h=h k +ofcp k and check the condition 
h, + a k p k eU, J(h, + a k p k )< J(h k ) (1) 
If this condition is fulfilled, then 

If this condition is not fulfilled we calculate the function J(h) at the point h=h k - 
<2kp k and check the condition 

. h k -a k p k eV, J(h k -a k p k )<J(h k ) (2) 
If this new condition is fulfilled, then 



N 



+ 1 
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Let us call a (k+l) xh iteration successful, if either condition (1) or (2) is fulfilled. If 
a (k+l) th iteration is not successful, then 

[ or* 5 otherwise 

Here A (0<X<1) is a parameter. The condition in (3) means that if in last N itera- 
5 tions with a step size parameter «k there is at least one successful iteration, then the step 
size parameter ak is kept for the following N iterations. If among the previous N itera- 
tions there is no successful one, then the step size parameter is reduced. 

In the described coordinate descent optimisation, if the function J(h) is convex and 
differentiate on U, then the sequence h k converges to the solution h for any initial val- 
10 ues ho in U, ao>0, and 0<A<L This result can be found, for example, in the book by F. 
P. Vasiliev, "Numerical methods for solving optimisation problems" (in Russian), 
Nauka, Moscow, 1988, page 345. 

In the LS problem the function J(h) is convex and differentiate; this is known, for 
example, from the book by F. P. Vasiliev, "Numerical methods for solving optimisation 
15 problems" (in Russian), Nauka, Moscow, 1988, page 169. Therefore, the coordinate de- 
scent optimisation when applied to the LS problem provides convergence to the solu- 
tion. The coordinate descent optimisation applied for finding the minimum of the func- 
tion J(h) in the LS problem is as follows. 

Let us pk=ei. At the (k+lf* iteration the conditions 
20 AJ(h k ) = J(h k + a k e, ) - J(h k ) < 0 

and 

AJ(h k ) = J{h k -a k e l )-J(h k )<0 

should be checked. Since only i' th element in the vector e, is not zero and the matrix 
R is symmetric, meaning that R(m,n)=R(n,m), we have 

25 AJ(h k ) = J(h k +a k e i )-J(h k ) = ^-\Q(i) + a k R(i,i)] 

where 

Q(i) = -2/?<7) + 2f> ( *>(m)tf(m, f) 

ro=l 

and h w (m). is m element of the vector hk. As ak>0 the condition 
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AJ(h k ) = J(h k +a k e ik )-J(h k )<0 

can now be represented as 
a k R(iJ)<-Q(i) 

Similarly, the condition 

AJ(h k ) = J(h k -a k e.)-J(h k )<0 

can be represented as 
a k R(U)>Q(i) 

Let Q k denote an auxiliary vector Q at k^ 1 iteration. If the (k+lf* iteration is not 
successful, then 

If the iteration is successful, then 

h (k+l) (i) = h (k) (i)±a k9 

Q {M \n) = Q w (n)±2a k R{nJ\ n = l 9 ... 9 N. 

where Q^(n) denotes n th element of the vector Qk. Initial values ho and Qo of the 
vectors h and Q can be chosen, for example, as follows: 

fc (0) («) = 0, n = l,...,N, 
Q (0 \n) = -2/3(n\ n = \^N. 

This means that the vector h is initialised by zeros, and the vector Q is initialised 
by the right side vector (5 of the system of linear equations Rh=p. 

If we choose H=2 q where q is an integer, ao=H/2, and A.=l/2, we avoid explicit 
multiplications by the step-size parameter a&. Such a choice allows multiplications to 
be implemented by means of bit-shift operations. Also, the described technique does not 
use any division operations. 

The method according to the invention may be characterised by the following 
parameters: (1) R is an NxN coefficient matrix of the system of linear equations; (2) p is 
a Nxl right-side vector of the system of linear equations; (2) h is a Nxl solution vector 
of the system of linear equations; (1) is the size of the solution vector; (2) M is the 
number of bits representing the solution; (3) Nit is the maximum number of passes for 
each bit of the solution; (4) H is a positive number defining the amplitude range of the 
solution. 



specification UNTVYORK-P0982.WOP-2002-04-06.doc 



- 13- 



For example, the method of solving systems of linear equations in accordance with 
the invention can be described as follows. 

Elements of the solution vector are represented as fixed-point binary words each 
consisting of at least one bit. An initialisation of the solution vector and an auxiliary 
vector are performed. The method performs bit-wise iterations for solving the system of 
linear equations, starting from the most significant bit and proceeding with the next less 
significant bit if a finishing condition is fulfilled. For each bit representing the binary 
words, the method performs at least one pass through all elements of the solution vec- 
tor, updating elements of the solution vector. The method stops solving the system of 
linear equations when a stopping condition is fulfilled. 

Elements of the solution vector are initialised as zeros; other initialisations are also 
possible. Elements of the auxiliary vector are initialised by corresponding elements of 
the right-side vector of the system of linear equations; other initialisations are also pos- 
sible. 

In the pass, in turn for each element of the solution vector a condition success- 
ful/unsuccessful is checked. If successful, the solution vector and the auxiliary vectors 
are updated. 

In the solution vector, the only element is updated, for which the condition success- 
ful/unsuccessful is checked. 

When checking the condition successful/unsuccessful for an element of the solu- 
tion vector, a finite number of possible updates of the element are checked to find a 
preferable update. The possible number of updates can be, for example, two (for real- 
valued data) or four (for complex-valued data); other variants are also possible. For 
each possible update an auxiliary value is calculated based on the corresponding ele- 
ment of the auxiliary vector and the minimum among the auxiliary values is found. The 
minimum is compared with a threshold calculated by using a diagonal element of the 
coefficient matrix, the diagonal element corresponding to the checked element of the 
solution vector, and a step-size parameter. The condition is successful if the minimum 
is less than the threshold, and it is unsuccessful if the minimum is higher than or equal 
to the threshold. The minimum indicates the preferable update and the updated element 
of the solution vector is replaced by the preferable update. For the successful condition, 
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the auxiliary vector is updated by using elements of a row of the coefficient matrix; this 
row corresponds to the updated element of the solution vector. 

When analysing an element of the solution vector the order of the analysis can be 
either arbitrary or a specific, for example, starting from an element whose position cor- 
responds to the position of an element of the auxiliary vector with a maximum ampli- 
tude and in the order of reducing the amplitude. In such a case, the method may also 
include an operation of ordering the elements of the auxiliary vector. 

The finishing condition is that in the last pass there was no successful update. 

The step-size parameter is decreased after analysing each bit of the solution; pref- 
erably, this decrease is by a factor of two. 

The method stops when a predefined number M of bit- wise iterations has been per- 
formed; this number defines valuable bits of the binary words representing elements of 
the solution vector. The method can stop if a condition that a predefined number Nit of 
passes through all elements of the solution vector is exceeded in a bit-wise iteration. 
Another condition to stop the method may be that the computer time intended for per- 
forming this method is finished. Other conditions may also be used. 

The solution can be obtained, for example, with a computer system programmed to 
perform a sequence of operations. 

1. Implementation of the method for real-valued data 

For real-valued data the method can be implemented by using the following se- 
quence of operations as represented in Figure 1 . 

Initialisation of the solution vector h and an auxiliary vector Q: h(r)=0, Q(r)=-j3(r) 3 
r=l, . . . , N Also, at this step, a bit number m, a pass number it, and an element number 
p are set to zero: m=0, it=0, p=0. 

Bit number is increased: m=m+l. 

A step-size parameter is calculated as d=2' m H. The preferable way to chose the 
amplitude range of the solution is H=2 q where q is an integer so that the step-size pa- 
rameter is d=2 q ' m and further operations involving d are performed without multiplica- 
tions. 

Pass number is increased: it=it+l. A binary parameter Flag is set to zero: Flag=0. 
Element number p is increased: p=p+L 
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An argument is found which indicates the minimum of the three values: Q(p), - 
Q(p) 9 and -(d/2)R(p 9 p): arg=min{Q(p), -Q(p), -(d/2)R(p,p)}. For calculating the third 
value, we avoid explicit multiplication because the step-size parameter d is a power of 
two and this multiplication is performed by a bit-shift operation. For finding the mini- 
5 mum we need to perform two operations of comparison, which require no multiplica- 
tion. 

The condition arg=],2, or 3 is checked. 

If arg=l then the following update of the element of the solution vector is per- 
formed: h(p)=h(p)+d. Also, all elements of the auxiliary vector Q are updated as 

1° Q( r ) = Q( r )+d-R(p,r) for r=l, ...,N; the binary parameter Flag is set to one: Flag=l, indi- 
cating that this iteration is successful. If arg=2 then the following update of the solution 
is performed: h(p)=h(p)-d. Also, all elements of the auxiliary vector Q are updated as 
Q( r )=Q( r )-d'R(P> r ) f° r r== l> — the binary parameter Flag is set to one: Flag=I, indi- 
cating that this iteration is successful. 

15 The condition p=N is checked. This equality means that in the current pass all N 

elements of the solution vector have been already analysed. If this condition is not satis- 
fied, the algorithm analyses the next element (step 5). 

If the condition p=N is satisfied the algorithm checks another condition: Flag=0. 
If the latter is not true, i.e., Flagyl meaning that in the course of the last pass, in- 

20 eluding sequential analysis of all TV elements of the solution vector, there is at least one 
successful update of elements of the solution vector h and the auxiliary vector Q, then a 
stopping condition "iterations? " is checked. This condition can be, for example, it=Nit, 
i.e., that the pass number it is equal to the maximum possible number of passes Nit 
which is a predefined parameter. This condition can also be based on checking the rest 

25 of a computer time predefined for solving the system of linear equations; other condi- 
tions can also be used. If this condition is satisfied then the method stops; in such a case 
the vector h contains N elements of the current solution. The bit number m characterises 
the accuracy of the solution. 

If the condition "iterations? " is not satisfied, a new pass is performed (step 4). 

30 If Flag=0, the condition m=M is checked; this condition means that all iWbits of 

the solution have been already analysed. If the condition m=M is not satisfied, the 
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method proceeds with analysing the next bit of the solution {step 2), i.e., performs a new 
bit-wise iteration. 

If m =M, the method stops; in such a case the vector h contains N elements of the 
solution. The predefined number of bits M characterises a predefined accuracy of the 
solution. 

2. Another implementation of the method for real-valued data 

For real-valued data the method can also be implemented by using an auxiliary 
delta array A and the following operations as represented in Figure 2. 

Initialisation of the solution vector h and an auxiliary vector Q: h(r)=0, Q(r)=-J3(r), 
r=l, . . . , K A bit number m, a pass number it, and an element number p are set to zero: 
m=0, it=0, p=0. 

Bit number is increased: m=m+l. 

A step-size parameter is calculated as d=2~ m H. The preferable way to chose the 
amplitude range of the solution is H=2 q where q is an integer so that the step-size pa- 
rameter is d=2 q ~ m and further operations involving d are performed without multiplica- 
tions. An delta array is initialised; this array has two elements: A(l)=d and A(2)~d. 

Pass number is increased: it=it+l. A binary parameter Flag is set to zero: Flag=0. 

Element number p is increased: p=p+l: 

An argument is found which indicates the minimum of the three values: Q(p), - 
Q(p)] and -(d/2)R(p i p): arg=min{Q(p), -Q(p), ~(d/2)R(p,p)}. For calculating the third 
value, we avoid explicit multiplication because the step-size parameter d is a power of 
two and this multiplication is performed by a bit-shift operation. For finding the mini- 
mum we need to perform two operations of comparison, which require no multiplica- 
tion. 

The condition arg< 3 is checked. 

If arg<3 then the following update of the solution is performed: h(p)^h(p)+ 
A(arg). Also, all elements of the auxiliary vector Q are updated as Q(r)=Q(r)+ 
A(arg) -R(p,r) for r=l, ...,N; the binary parameter Flag is set to one: Flagyl, indicating 
that this iteration is successful. 

The condition p=N is checked. This equality means that in the current pass all N 
elements of the solution vector have been already analysed. If this condition is not satis- 
fied, the next element of the solution vector is analysed {step 5). 



specification UNFVYORK~P0982.WOP-2002-04-06.doc 



- 17- 



If the condition p=N is satisfied, another condition is checked: Flag=0. 
If the latter is not true, i.e., Flag=], meaning that in the course of the last pass, in- 
cluding sequential analysis of all N elements of the solution vector, there is at least one 
update of elements of the solution vector h and the auxiliary vector Q, then a stopping 
5 condition "iterations?" is checked. This condition can be, for example, it=Nit, i.e., that 
the pass number it is equal to the maximum possible number of passes Nit which is a 
predefined parameter. This condition can also be based on checking the rest of a com- 
puter time intended for solving the system of linear equations; other conditions can also 
be used. If this condition is satisfied the methods stops; in such a case the vector h con- 
10 tains N elements of the solution. The bit number m characterises the accuracy of the so- 
lution. 

If the condition "iterations? " is not satisfied, a new pass is performed {step 4). 
If Flag=0, the condition m=M is checked; this condition means that all Mbits of 
the solution have been already analysed. If the condition m=M is not satisfied, a new 
15 bit- wise iteration is performed by analysing the next bit of the solution (step 2). 

If m=M, the methods stops; in such a case the vector h contains N elements of the 
solution vector with the predefined number M of valuable bits. The predefined number 
of bits M characterises a predefined accuracy of the solution. 
3. Program implementing the method for real-value data 
20 An example of implementation of the method in accordance with the invention for 

real-valuable data with an auxiliary delta array by a MATLAB program is as follows: 

function h=LS_real (R, N, M, H, Nit , beta ) ; 
h=zeros (N, 1 ) ; 
Q=-beta; 
25 d=H/2; 

for m=l:M 

delta (l)=d; delta (2)=-d; 

d=d/2; 

for it=l:Nit 
30 Flag=0; 

for p=l:N 

[val arg] =min ( [Q(p) ,-Q(p) , -abs (R(p,p)*d) ] ) ; 
if arg<3 
Flagyl; 

35 h (p) =h (p) +delta (arg) ; 
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Q=Q+delta(arg)*R(p, :); 

end 

end 

if Flag==0 break; end 

end 

end 

4. Another implementation of the method for real-valued data 

In some cases it is preferable to obtain such a solution of the system of the equa- 
tions that the amplitude range for the solution is predefined, i.e., there are lower hi and 
upper h2 bounds for elements of the solution vector; for example, hj=0 which means 
that the elements of the solution vector should be positive. For real-valued data the 
method can be implemented by using the following operations as represented in Figure 
3. 

Initialisation of the solution vector h and an auxiliary vector Q: h(r)=0, Q(r)=-/3(r), 
r=l, . . . , K A bit number m, a pass number zY, and an element number p are set to zero: 
m=0, it=0, p=0. 

Bit number is increased: m=m+l. 

A step-size parameter is calculated as d=2' m H. The preferable way to chose the 
amplitude range of the solution is H=2 q where q is an integer so that the step-size pa- 
rameter is d=2 q ' m and further operations involving d are performed without multiplica- 
tions. Ah delta array is initialised; this anSy has two elements: A'(l)=d and A(2)~d. 

Pass number is increased: it=it+l. A binary parameter Flag is set to zero: Flag=0. 

Element number p is increased: p=p+J. 

An argument is found which indicates the minimum of the three values: Q(p), - 
Q(p), and -(d/2)R(p i p): arg=min{Q(p), -Q(p), -(d/2)R(p,p)}, For calculating the third 
value, we avoid explicit multiplication because the step-size parameter d is a power of 
two and this multiplication is performed by a bit-shift operation. For finding the mini- 
mum we need to perform two operations of comparison, which require no multiplica- 
tion. 

The condition arg< 3 is checked. 

If arg<3 then the following update of the solution is performed: h(p)=h(p) + 
A(arg). Then the condition hj^i(p)<h 2 is checked. If it is true, then all elements of the 
auxiliary vector Q are updated as Q(r)=Q(r)+ A(arg) -R(p,r) for r=l,...,N; the binary 
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parameter Flag is set to one: Flag=l, indicating that this iteration is successful. If the 
condition hj<h(p)<h2 is not satisfied then the following update of the solution is per- 
formed: h(p)=h(p)- A(arg). Note that a more general approach can be implemented at 
this step when the condition h/(p)<h(p)<h2(p) is checked; here, hi(p) and h 2 (p) are ele- 
ments of two predefined Nxl vectors containing lower and upper bounds which are in- 
dividual for each element of the solution vector. 

The condition p=N is checked. This equality means that in the current pass all TV 
elements of the solution vector have been already analysed. If this condition is not satis- 
fied, the next element is analysed (step 5). 

If the condition p=N is satisfied, then another condition is checked: Flag=0. 

If the latter is not true, i.e., Flag=], meaning that in the course of the last iteration, 
including sequential analysis of all TV elements of the solution, there is at least one suc- 
cessful update of elements of the solution vector h and the auxiliary vector Q, then a 
stopping condition "iterations? " is checked. This condition can be, for example, it=Nit, 
i.e., that the pass number it is equal to the maximum possible number of passes Nit 
which is a predefined parameter. This condition can also be based on checking the rest 
of a computer time intended for solving the system of linear equations; other conditions 
can also be used. If this condition is satisfied the method stops; in such a case the vector 
h contains TV elements of the solution. The bit number m characterises an accuracy of 
the solution. 

if the condition "iterations? " is not satisfied, a new pass is performed (step 4). 
If Flag=0, the condition m=M is checked; this condition means that all Mbits of 
the solution have been already analysed. If the condition m=M is not satisfied, a new 
bit-wise iteration is performed by analysing the next bit of the solution (step 2). 

If m=M, the method stops; in such a case the vector h contains TV elements of the 
solution. The predefined number of bits M characterises a predefined accuracy of the 
solution. 

5. A general implementation of the method for complex-valued data 

In order to solve a linear system of equations Rh=P with complex-valued data one 
can use implementations of the method for real-valued data (for example, that described 
above) in application to a new real-valued system Ab=c with the following well-known 
replacements: the size of the real- valued system of linear equations is 2Nx2N (A is a 
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2Nx2N real -valued coefficient matrix, b is a 2Nxl real-valued solution vector, and c is a 
2Nxl real-valued right-side vector); elements of the vector b are defined as b(2n- 
l)=Re{h(n)} and b(2n)=Im{h(n)}\ elements of the vector c are c(2n-l)=Re{/3(n)} 
c(2n)=Im{/3(n)}\ elements of the matrix A are A(2m-l,2n-l)=Re{R(m,n)}, 
5 A(2m,2n)=Re{R(m,n)}, A(2m-l,2n)=Im{R(m,n)}, and A(2m,2n-l)=-Im{R(m,n)}; Re{} 
and Im{} denote, respectively, real and imaginary parts of a complex number. 
6. Implementation of the method for complex-valued data 

For complex-valued data the method can also be implemented by using the follow- 
ing sequence of operations as represented in Figure 4. 
10 Initialisation of the solution vector h and an auxiliary vector Q: h(r)=0, Q(r)=-/3(r), 

r=l, . . . , K A bit number m, a pass number it, and an element number p are set to zero: 
m=0, it=0, p=0. 

Bit number is increased: m=m+l. 

A step-size parameter is calculated as d=2' m H, The preferable way to chose the 
15 amplitude range of the solution is H=2 q where q is an integer so that the step-size pa- 
rameter is d=2 g ' m and further operations involving d are performed without multiplica- 
tions. 

Pass number is increased: it=it+l. A binary parameter Flag is set to zero: Flag=0. 
Element number p is increased: p=p+L 
20 An argument is found which indicates the minimum of the five values: Re[Q(pJ), - 

Re[Q<p)L MQ(p]h -Im[Q(p)L and ~(d/2)R(p,p): arg=min{Re[Q(p]), -Re[Q(p]), 
Im[Q(p]), -Im[Q(p)], -(d/2)R(p,p)}. For calculating the last value, we avoid explicit 
multiplication because the step-size parameter d is a power of two and this multiplica- 
tion is performed by a bit- shift operation. For finding the minimum we need to perform 
25 operations of comparison, which require no multiplication. 
The condition arg=l, 2,3,4 or 5 is checked. 

If arg=l then the following update of the solution is performed: h(p)=h(p)+d\ the 
latter means that the real part of the element h(p) is increased by d. Also, all elements of 
the auxiliary vector Q are updated as Q(r)=Q(r)+d-R(p,r) for r=7, ...,N; the binary pa- 
30 rameter Flag is set to one: Flag=l, indicating that the current pass is successful. If 
arg=2 then the following update of the solution is performed: h(p)=h(p)-d; the latter 
means that the real part of the element h(p) is decreased by d. Also, all elements of the 
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auxiliary vector Q are updated as Q(r)=Q(r)-d-R(p,r) for r=l, ...,N; the binary parameter 
Flag is set to one: Flagyl, indicating that the current pass is successful. If arg=3 then 
the following update of the solution is performed: h(p)=h(p)+i-d where i 2 ~l\ the latter 
means that the imaginary part of the element h(p) is increased by d. Also, all elements 
5 of the auxiliary vector Q are updated as Q(r)=Q(r)+i-dR(p,r) for r=l,...,N; the binary 
parameter Flag is set to one: Flag=l, indicating that the current pass is successful. If 
arg=4 then the following update of the solution is performed: h(p)=h(p)-i-d\ the latter 
means that the imaginary part of the element h(p) is decreased by d. Also, all elements 
of the auxiliary vector Q are updated as Q(r)=Q(r)-i-d'R(p,r) for r=l,...,N; the binary 
10 parameter Flag is set to one: Flag=l, indicating that the current pass is successful. 

The condition p=N is checked. This equality means that in the current pass all N 
elements of the solution vector have been already analysed. If this condition is not satis- 
fied, the next element is analysed (step 5). 

If the condition p=N is satisfied, another condition is checked: Flag=0. 
15 If the latter is not true, i.e., Flag=l, meaning that in the course of the last pass, in- 

cluding sequential analysis of all N elements of the solution, there is at least one suc- 
cessful update of elements of the solution vector h and the auxiliary vector Q, then a 
stopping condition "iterations?" is checked. This condition can be, for example, it=Nit, 
i.e., that the pass number it is equal to the maximum possible number of passes Nit 
20 which is a predefined parameter. This condition can also be based on checking the rest 
of a computer time intended for this algorithm; other conditions can also be used. If the 
stopping condition is satisfied the method stops; in such a case the vector h contains N 
elements of the current solution. The bit number m characterises the accuracy of the so- 
lution. 

25 If the condition "iterations? " is not satisfied, a new pass is performed (step 4). 

If Flag=0, the condition m=Mis checked; this condition means that all JWbits of 
the solution have been already analysed. If the condition m=M \s not satisfied, the next 
bit-wise iteration is performed by analysing the next bit of the solution (step 2). 

If the method stops; in such a case the vector h contains N elements of the 

30 current solution. The predefined number of bits M characterises an accuracy of the solu- 
tion. 
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In step 6 of the implementation described above, finding arg=min{Re[Q(p]), - 
Re[Q(p])> I™[Q(P])> -I m tQ(P)]> -(<V 2 )R(P>P)} is required. This can be implemented by 
using the following operations as represented in Figure 5. 

Input data are received: a=Re[Q(p)] and b=Im[(Q(p)] . 
5 The processing of the input data a and b is performed in parallel as follows. The 

condition a>0 is checked. If this condition is satisfied, then a=-a and a binary parame- 
ter I a =l. If the condition a>0 is not satisfied, then I a =0. The condition b>0 is checked. 
If this condition is satisfied, then b =-b and a binary parameter I b =L If the condition 
b>0 is not satisfied, then 
10 The condition a>b is checked. If it is satisfied, then h=l 9 and c=b. If the 

condition is not satisfied then Ii=I a h = 0, and c=a. 

The condition c>-(d/2)-R(p,p) is checked. If it is true, then I 3 =J; if it is false, then 
h=0. 

The binary values //, 1% h indicate the position of the minimum: I 3 =0/J means that 
15 (an update)/(no update) of the solution is required; I 2 =0/1 means that the real/imaginary 
part of the solution should be updated; Ij=0/I means that the update should be posi- 
tive/negative. 

7. Another implementation of the method for complex-valued data 

For complex-valued data the method can also be implemented by using an auxil- 
20 iary delta array A and the following operations as represented in Figure 6. 

Initialisation of the solution vector h and an auxiliary vector Q: h(r)=0 } Q(r)=-(3(r), 
r=7, . . . , N. A bit number m, a pass number it, and an element number p are set to zero: 
m=0, it=0, p=0. 

Bit number is increased: m=m+l. 
25 A step-size parameter is calculated as d=2' m H. The preferable way to chose the 

amplitude range of the solution is H=2 q where q is an integer so that the step-size pa- 
rameter is d=2 q ' m and further operations involving d are performed without multiplica- 
tions. An delta array is initialised; this array has four elements: A(l)=d } A(2)~d, 
A(3)=i-d, andA(4)=-i~d. 
30 Pass number is increased: it=it+l. A binary parameter Flag is set to zero: Flag=0. 

Element number p is increased: p=p+l* 
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An argument is found which indicates the minimum of the five values: Re[Q(p])> - 
RefQfr)], Im[Q(p])> -Im[Q(p)h and -(d/2)R(p,p): arg=min{Re[Q(p]), -RefQfr]), 
I m [Q(p])> -I m [Q(p)]> -(d/2)R(p,p)}. For calculating the last value, we avoid explicit 
multiplication because the step-size parameter d is a power of two and this multiplica- 
tion is performed by a bit-shift operation. For finding the minimum we need to perform 
operations of comparison, which require no multiplication. 

The condition arg< 5 is checked. 

If arg<5 then the following update of the solution is performed: h(jp)=h(p) + 
A(arg). Also, all elements of the auxiliary vector Q are updated as Q(r)=Q(r)+ 
A(arg) -R(p,r) for r=l, ...,N; the binary parameter Flag is set to one: Flag=l, indicating 
that the current iteration is successful. 

The condition p=N is checked. This equality means that in the current pass all N 
elements of the solution vector have been already analysed. If this condition is not satis- 
fied, the next element is analysed (step 5). 

If the condition p=N is satisfied, another condition is checked: Flag=0. 

If the latter is not true, i.e., Flag=l, meaning that in the course of the last pass, in- 
cluding sequential analysis of all N elements of the solution, there is at least one suc- 
cessful update of elements of the solution vector h and the auxiliary vector Q, then a 
stopping condition "iterations? " is checked. This condition can be, for example, it=Nit, 
i.e., that the pass number it is equal to the maximum possible number of passes Nit 
which is a predefined parameter. This condition can also be based on checking the rest 
of a computer time intended for this algorithm; other conditions can also be used. If this 
condition is satisfied the method stops; in such a case the vector h contains N elements 
of the current solution. The bit number m characterises an accuracy of the solution. 

If the condition "iterations? " is not satisfied, a new pass is performed (step 4). 

If Flag=0, the condition m=M is checked; this condition means that all M bits of 
the solution have been already analysed. If the condition m=M is not satisfied, the next 
bit-wise iteration is performed by analysing the next bit of the solution (step 2). 

If m=M 9 the method stops; in such a case the vector h contains N elements of the 
solution. The predefined number of bits M characterises an accuracy of the solution. 
8. Program implementing the method for complex-valued data 

An example of implementation of the described method for complex-valued data 
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with an auxiliary delta array by a MATLAB program is as follows: 

function h=LS_cmplx (R, N, M, H, Nit, beta) ; 
h=zeros (N, 1 ) ; 
Q=-beta; 
5 • d=H/2; 

for m=l:M 

delta (l)=d; delta(2)=-d; delta ( 3 ) =i*d; delta ( 4 ) =-i*d; 
d-d/2; 

for it=l:Nit 
10 Flag=0; 

for p=l:N 

[val arg]=min( [real (Q(p) ), -real (Q(p) ) , . . - 

imag{Q(p) ) f -imag(Q(p) ) ,-R(p,p) *d] ) ; 

if arg<5 
15 Flagyl; 

h (p) =h (p) +delta (arg) ; 
Q=Q+delta (arg) *R(p, : ) ; 

end 

end 

20 if Flag==0 break; end 

end 

end 

Figure 7 shows a computer system for solving systems of linear equations by the 
method according to the invention. The computer system comprises a parameter inter- 
25 face and a computer processor. The parameter interface receives parameter signals rep- 
resenting the coefficient matrix R and the vector P of the system of linear equations 
Rh=p. The output of the computer processor comprises elements of the solution vector 
h. 

9. Computer system 

30 Figure 8 represents a computer system in accordance with the invention, the com- 

puter system implementing the method of solving systems of linear equations in accor- 
dance with the invention. The computer system comprises a Host processor, a Block of 
h-updates, a Block of Q-updates, a Block of R-memory, a Block of minimisation, and a 
Controller; it also contains a Host Bus and an Internal Bus. The Host processor is con- 

35 nected with the Block of h-updates, the Block of Q-updates, the Block of R-memory, 
and the Controller by means of the Host Bus. The Controller is connected with the 
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Block of h-updates, the Block of Q-updates, the Block of R-memory, and the Block of 
minimisation by means of the Internal Bus. In addition, the Block of Q-updates is con- 
nected with the Block of minimisation, and the Block of R-memory is connected with 
the Block of Q-updates and the Block of minimisation. 

The Host processor produces itself or receives from other devices (not shown here) 
parameter signals representing the coefficient matrix R and the vector p of the linear 
system Rh=p. 

There are three stages of operation of the computer system: an initialisation stage, a 
calculating stage, and a stage of reading results. 

At the initialisation stage, the Host processor first sends (through the Host Bus) to 
all blocks a signal Init meaning the beginning of the initialisation stage, and then the 
Host processor sends signals representing the vector p to the Block of Q-updates and 
signals representing the coefficient matrix R to the Block of R-memory. In some im- 
plementations, the signals representing the coefficient matrix R can be stored in the 
Block of R-memory before the initialisation stage. Also, the Host processor sends sig- 
nals representing parameters of the method through the Host Bus to the other blocks, 
preferably to the Controller; these parameters may be the size of the solution vector (N), 
number of bits representing the solution (M), maximum number of passes (Nit), ampli- 
tude range of the solution (H). Also, among these parameters may be a maximum opera- 
tion time, i.e., a time duration during which the calculation of the solution should be 
done. In some implementations, these parameters can be stored in the Controller instead 
of sending them by the Host processor. Upon receiving the signal Init the Block of h- 
updates sets elements of the vector h to zero. 

At the calculating stage, the Controller forms control signals needed for performing 
the method and sends them to the other blocks through the Internal Bus and to the Host 
processor through the Host Bus. The control signals may comprise the bit number m, 
the pass number it, the element number p, a number r of the element to be updated, sig- 
nals representing commands to begin and stop operations, address signals to memory 
blocks; other control signals can also be sent by the Controller to the blocks. The Block 
of h-updates updates the solution vector h. The Block of Q-updates updates the auxil- 
iary vector Q. The Block of R-memory sends signals representing the coefficient matrix 
R to the Block of Q-updates and the Block of minimisation according to address signals 
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received from the Controller. The Block of minimisation receives (from the Block of Q- 
updates) signals representing elements of the auxiliary vector Q and (from the Block of 
R-memory) elements of the matrix R and finds a position (arg) of the minimum re- 
quired for controlling the updates of the vectors h and Q. Signals representing the pa- 
5 rameter arg are sent to the Internal Bus, from where they are received by the Controller 
and may also be received by the other blocks. 

At the stage of reading results the Host processor forms address signals for the so- 
lution vector h 5 sends these signals to the Block of h-updates and reads signals repre- 
senting the solution vector h. 
10 In some implementations, the Host processor can stop solving the linear system at 

any moment and read the current solution from the Block of h-updates. It can also read 
other parameters, such, for example, as the current bit number m. The latter allows the 
Host processor to obtain information of the accuracy of the current solution. 

In some implementations, functions of the Controller may be performed by the 
15 Host processor; in such a case the Host Bus and the Internal Bus may be connected. In 
some implementations, the Block of minimisation may be included in the Block of Q- 
updates. Other reconfigurations of the computer system are also possible. 

A variant of implementation of the block of R-memory is shown in Figure 9. 

At the initialisation stage, the block receives (from the Host Bus) the signal Init, 
20 "indicating the initialisation stage, and signals representing row (Init row address) and 
column (Init column address) numbers of elements of the coefficient matrix R as well 
as elements of the matrix {Coefficients of the matrix R). The signal Init at the address 
inputs of the Row and Column multiplexers allows the multiplexers to transmit the ad- 
dress signals to the address input of a Memory block of the coefficient matrix R. The 
25 coefficients are recorded to the Memory block of the coefficient matrix R. In another 
implementation, the Memory block of the coefficient matrix R can be a read-only- 
memory (ROM) with a predefined matrix R; in such a case, the Row and Column mul- 
tiplexers should be excluded (together with corresponding signals from the Host Bus). 

At the calculating stage, the block receives (from the Internal Bus) signals repre- 
30 senting an element number p, a number r of the element to be updated, the step-size pa- 
rameter d, and a control signal Update/Analysis; other signals to control (reading 
from)/(writing to) the Memory block of the coefficient matrix R can also be received 
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from the Internal Bus. There are two operation modes of the block in the calculating 
stage: Update and Analysis modes. When in the Analysis mode, the row and column ad- 
dresses of the Memory block of the coefficient matrix R are equal to the element num- 
ber p., i.e., p th diagonal element of the matrix R is read from the Memory block of the 
coefficient matrix R. This is provided by the Update multiplexer whose address input is 
connected with the signal Update/Analysis. When in the Update mode, the row address 
of the Memory block of the coefficient matrix R is r. Output of the Memory block of 
the coefficient matrix R is bit-shifted according to the step-size parameter d by means of 
a Bit shift block. In the Update mode, output data of the block are sent to the Block of 
Q-updates; in the Analysis mode, output data of the block are sent to the Block of mini- 
misation. 

At all these stages the block can also receive (from the Host Bus and/or from the 
Internal Bus) other signals to control the Memory block of the coefficient matrix R. 

A variant of implementation of the block of h-updates is shown in Figure 10; this 
variant is for complex-valued data. 

At the initialisation stage, the block receives the signal Init from the Host Bus; 
upon this signal all elements of a Memory block of the solution vector h are set to zero. 

At the calculating stage, the block receives (from the Internal Bus) signals repre- 
senting the element number p, the step-size parameter d, and signals // and I 2 indicating 
whether the real or imaginary part of the element h(p) is updated (I2) and whether addi- 
tion or subtraction of the step-size parameter d is performed (//). Also, a signal Up- 
date/Reading defines that an Update/Reading multiplexer transmits signals /? and p to 
the address input of the Memory block of the solution vector h. The real or imaginary 
part (depending of the signal I 2 ) of the parameter h(p) is read from the Memory block of 
the solution vector h. In an Adder-subtractor, this part of the parameter h(p) is updated 
by adding or subtracting (depending on the signal I } ) the step-size parameter d. The up- 
dated parameter h(p) is re-written in the Memory block of the solution vector h at the 
same address. 

At the stage of reading the results, the block receives (from the Host Bus and/or 
may be from the Internal Bus) signals representing an address of reading from the 
Memory block of the solution vector h (Reading address) and the signal Up- 
date/Reading. The signal Reading address is passed to the address input of the Memory 
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block of the solution vector h, and output signals {Elements of the solution vector) rep- 
resenting the solution vector h are sent to the Host Bus. 

At all these stages the block can also receive (from the Host Bus and/or from the 
Internal Bus) other signals to control the Memory block of the solution vector h. 
5 A variant of implementation of the Block of Q-updates is shown in Figure 1 1 ; this 

variant is for complex-valued data. 

At the initialisation stage, the block receives (from the Host Bus and/or from the 
Internal Bus) a signal Update/Init, which allows the initialisation data (signals repre- 
senting elements of the vector (3) to pass through a Data multiplexer and to be written in 

10 the Memory block of the auxiliary vector Q. Also, the block receives (from the Host 
Bus) signals Init address representing an address and a signal Update/ Analysis/Init al- 
lowing the address Init address to pass through an Update/Analysis/Init multiplexer to 
the address input of the Memory block of the auxiliary vector Q. 

At the calculating stage, the block receives (from the Internal Bus) signals repre- 

15 senting the element number p, the number r of the element to be updated, the step-size 
parameter d, and signals // and h indicating a type of the update of the vector Q. Also, 
the block receives (from the Block of R-memory) signals representing the value 
d-R(p,r). At this stage, the block operates in one of the following modes: Update mode 
or Analysis mode. In the Analysis mode, the Update/Analysis/Init multiplexer passes 

20 signals representing the element number p to the address input of the Memory block of 
the auxiliary vector Q; the Memory block of the auxiliary vector Q sends signals repre- 
senting the element Q(p) to the Block of minimisation. In the Update mode, the Up- 
date/Analysis/Init multiplexer passes signals representing the number r of the element to 
be updated to the address input of the Memory block of the auxiliary vector Q; the 

25 Memory block of the auxiliary vector Q sends signals representing the element Q(r) to a 
Complex adder-subtractor. Simultaneously, signals representing the value d-R(p y r) are at 
another data input of the Complex adder-subtractor. The signals Ij and h define the type 
of update: if 1 ) =0/1 the Complex adder-subtractor performs summation/subtraction; the 
signal h=0/l indicates that the value d-R(p,r) is (not multiplied)/(multiplied) by i; the 

30 latter means that real and imaginary parts of the value are replaced by each other. The 
updated element Q(r) at the output of the complex Adder-subtractor is re- written to the 
Memory block of the auxiliary vector Q at the same address. 
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At all these stages the block can also receive (from the Host Bus and/or from the 
Internal Bus) signals to control the Memory block of the auxiliary vector Q. 

A variant of implementation of the Block of minimisation is shown in Figure 12; 
this variant is for complex-valued data. The block receives (from the Block of Q- 
updates) real and imaginary parts of a p th element of the auxiliary vector Q: a=Re[Q(p)] 
and b=Im[Q(p)]. The block also receives (from the Block of R-memory) a signal repre- 
senting the parameter d>R(p,p). A Converter-a converts the signal a to two signals repre- 
senting, respectively, the sign sgn(a) and module \a\ of a. Similarly, a Converter-b con- 
verts b to sgn(b) and \b\. .First comparator compares \a\ and \b\ and forms a binary signal 
h=0/l if |a|>|A| or |a|<=|6| 5 respectively. A multiplexer passes the maximum of \a\ and 
| b | to a Second comparator. At another input of the Second comparator are signals rep- 
resenting the parameter d-R(p,p). The second comparator compares the input signals and 
forms a binary signal I 3 =0/1 if max(\a\,\b\)> d-R(p,p) or max(\a\,\b\)<= d-R(p,p), re- 
spectively. The signals representing sgn(a) and sgn(b) go to the data inputs of a Bit mul- 
tiplexer whose address input is connected to the signal h\ output of the Bit multiplexer 
forms a signal //. 

The equality h=0 means that the current iteration is successful and the element 
h(p) of the solution vector h should be updated as well as all elements of the auxiliary 
vector Q should be updated, while 7?=7 means that the current iteration is not successful 
and no update is required. The equality I 2 =0/1 means that the real/imaginary part of the 
element h(p) of the solution vector h should be updated; it also indicates that the value 
d-R(p,r) is (not multiplied)/(multiplied) by i\ the latter means that real and imaginary 
parts of the value are replaced by each other. The equality Ij =0/1 means that addi- 
tion/subtraction is used when updating. 

The signals Ij, h and Is are sent to the Internal Bus. From the Internal Bus, these 
signals are received by the Controller, Block of h-update, and Block of Q-updates. In 
some implementations, these signals can be modified by the Controller before the use in 
the Block of h-updates and Block of Q-update. 
10. Multiuser receiving method and multiuser receiver 

The invention is also applicable to all kinds of CDMA data transmission systems. 
In the following, the invention will be described in greater detail in a cellular radio sys- 
tem without, however, limiting it thereto. In a cellular radio system, in which the 

specification UNIVYORK~P0982. WOP-2002-04-06.doc 



-30- 



method of multiuser receiving according to the invention can be applied, each cell com- 
prises at least one base station communicating with mobile stations. The base station 
transmits calls from the mobile stations to a public telephone network or to another mo- 
bile station. A succession of data signals of a mobile user, translated into bits, is trans- 
5 mitted at a rate of a plurality of chips per bit, spread by a respective spreading code. All 
the mobile stations transmit at the same frequency to the base station, which distin- 
guishes the transmissions of different mobile stations on the basis of respective different 
spreading codes. The signals of the mobile stations can also propagate along several dif- 
ferent paths to the base station. 

10 In the multiuser receiving method in accordance with the invention, a signal re- 

ceived by an antenna in a receiver is transmitted through radio frequency parts and other 
necessary signal processing parts, to filters matched with the spreading codes used, and 
in the filters the signal can be restored to the original band. The matched filters output 
the signals of the respective users with the respective spreading codes extracted. The 

15 signals at the matched filter outputs are grouped in a Nxl vector p, where N is the num- 
ber of the spreading codes. Results of the matched filtering are transformed to a solution 
vector h whose values are represented by a solution of a system of linear equations of 
the kind Rh=P where R is a NxN cross-correlation matrix of the spreading codes and h 
is the Nxl solution vector. When transforming the matched-filter outputs, the method of 

20 " solving systems of linear equations according to the invention is applied. The method of 
solving systems of linear equations can be implemented as a sequence of operations in a 
computer system. In another embodiment of the invention, when transforming the 
matched-filter output vector, the computer system according to the invention is applied. 
The signals representing elements of the solution vector h are then forwarded to a detec- 

25 tor, where a data decision is made. The cross-correlation matrix of the spreading codes 
comprises computed cross-correlations of the spreading codes used, and it indicates 
how much the signals of different users interfere with one another. Since the solving of 
the system of linear equations is performed by the method or the computer system ac- 
cording to the invention, which can be implemented without multiplication and division 

30 operations, the multiuser receiving is simplified. 

The multiuser receiver in accordance with the invention can be implemented as 
represented in Figure 13. Figure 13 is a block diagram illustrating the basic structure of 
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the multiuser receiver according to the invention. The multiuser receiver comprises an 
antenna with which a received signal is forwarded to radio frequency (RF) parts, where 
the radio frequency signal is converted to an intermediate frequency. From the radio 
frequency parts, the signal is supplied to an analog-to-digital converter (A/D), where the 
received analog signal is converted to digital form. The converted digital signal is sup- 
plied to detector means (DET), and the output signal is then forwarded to the other parts 
of the receiver. The receiver also comprises control means (CNTL), which control the 
operation of the above-mentioned blocks. 

Figure 14 illustrates in greater detail implementation of the detector means accord- 
ing to one preferred embodiment of the invention. The detector block comprises filters 
matched with the spreading codes, in the matched filters, the spread spectrum signal is 
restored to the original narrow band. The number of matched filters is equal to the num- 
ber of signals that are to be used in the detection, such as the number of users N. Signals 
obtained from the matched filters are supplied to a computer system and transformed 
into a solution vector h whose elements are represented by a solution of a system of lin- 
ear equations of the kind Rh=p where R is a NxN cross-correlation matrix of the spread- 
ing codes and h is the Nxl solution vector. When transforming the matched-filter output 
vector, the method of solving systems of linear equations according to the invention is 
applied; the method of solving systems of linear equations can be implemented as a se- 
quence of operations in a computer system. In another embodiment of the invention, 
when transforming the signals obtained from the matched-filter, the computer system 
according the invention can be applied. The signals representing elements of the vector 
h are then forwarded to a means making a data symbol decision. The receiver naturally 
also comprises other components, such as filters or means for calculating the cross- 
correlation matrix, but for the sake of clarity and since they are not essential to the pre- 
sent invention, they have not been shown. 
11. Adaptive filtering 

As an example, we consider adaptive filtering in application to echo cancellation in 
a transmission system. However, other applications of adaptive filtering according to 
the invention are also possible. Prior to description of the preferred embodiment, de- 
scription will be made as to general arrangement and operation of the data transmission 
system with an echo canceller so as to help better understanding of he present invention. 
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Referring to Figure 15 5 an incoming signal propagates as a first signal x(t) through 
a first transmission line to an unknown system which produces a second signal y(t) 
transmitted to a second transmission line. Ideally, the second signal y(t) contains only a 
signal of interest s(t)\ y(t)=s(t), and does not contain any echo: e(t)=0. In practice, how- 
ever, a part of the first signal x(t) leaks to the second transmission line through the un- 
known system. Then, the signal of interest s(t) is mixed with the echo to form the sec- 
ond signal: y(t)=s(t)+e(t). The echo eft) should be cancelled from the second signal so 
as to transmit the second signal as the signal of interest s(t) alone. The unknown system 
can be a hybrid transformer; a transmission system with a hybrid transformer and echo 
cancellation is discussed in greater detail in the Patent US 5062102, "Echo canceller 
with means for determining filter coefficients from autocorrelation and cross-correlation 
coefficients" by Tetsu Taguchi. The unknown system can also be an acoustic system as, 
for example, is discussed in greater detail in Patent WO 00/38319, "Stable adaptive fil- 
ter and method" by Ding Heping. 

In order to cancel the echo from the second signal, the echo canceller comprises a 
filter such as a transversal filter coupled to the first transmission line and responsive to 
the first signal x(t) for producing an echo replica or an estimated echo signal e\(t) de- 
termined in accordance with the filter coefficients h(i) (1=1,2, ...,N). A subtracter is cou- 
pled to the filter and connected in the second transmission line for subtracting the esti- 
" mated echo signal e } (t) from the second signal y(t)=s(t)+e(t) on the second transmission 
line so that the echo is reduced or cancelled: yi(t)=y(t)-e I (t)=s(t)+[e(t)-e 1 (t)] \ The echo 
canceller also comprises a coefficient generating circuit for generating the filter coeffi- 
cients h(i) so as to determine the estimated echo signal ej(t). If the coefficient generat- 
ing circuit can generate the optimum filter coefficients so as to provide e(t)=ej(t), the 
echo is completely cancelled. When e(t)*ej(t) 9 [e(t)-ej(t)J still remains as a residual 
echo. That is, the echo cancellation is defined by performance of the coefficient generat- 
ing circuit. 

In the shown embodiment in Figure 16, the coefficient generating circuit comprises 
a first, a second, and a third circuit. The first circuit is coupled to the first transmission 
line and calculates autocorrelation coefficients of the first signal x(t). The second circuit 
is coupled to the first and the second transmission lines and calculates cross-correlation 
coefficients between the first signal x(t) and the second signal y(t). Accordingly, the first 
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and the second circuits are an autocorrelation coefficient calculator and a cross- 
correlation coefficient calculator, respectively. 

The third circuit is coupled to the autocorrelation coefficient calculator and the 
cross-correlation coefficient calculator and decides the filter coefficients h(i) 
(i=l,2, ... f N). The third circuit is called a filter coefficient decision circuit. It is generally 
known in the prior art that the optimal filter coefficients are that obtained by solving a 
system of linear equations of the kind Rh=p where R is a NxN coefficient matrix com- 
prising the autocorrelation coefficients R(m y n) of the first signal x(t), 

T-\n-m\ 

R{m,ri)= ^x{t)x(t+\m-n\\ m,n = 1,2,..., TV, 

where /=/, are discrete time moments, a Nxl vector p comprises the cross- 
correlation coefficients, 

/=i 

and h is a Nxl solution vector comprising the filter coefficients h(i) (i=l,2, ...,N)l 
this is discussed in greater detail in the Patent US 5062102, "Echo canceller with means 
for determining filter coefficients from autocorrelation and cross-correlation coeffi- 
cients" by Tetsu Taguchi. 

The autocorrelation coefficients and the cross-correlation coefficients are provided 
to the third" circuit from the autocorrelation coefficient calculator and the cross- 
correlation coefficient calculator, respectively. The third circuit in the embodiment is a 
computer system for solving systems of linear equations. When the computer system is 
implemented by using the described above method for solving systems of linear equa- 
tions in accordance with the invention, the convergence speed increases and truncation 
errors are eliminated. As a result, the echo canceller exploiting such a computer system 
allows better cancellation of the echo. In another embodiment, the computer system is 
implemented as the computer system in accordance with the invention; this allows the 
filter coefficients to be calculated without multiplication or division operations, leading 
to simpler implementation. The latter allows the sampling frequency of the adaptive fil- 
ter to be increased and to be applied in data transmission systems with higher data rates. 
This also allows the number of filter coefficients N to be increased, providing a better 
echo cancellation, i.e., a smaller residual echo. 
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Although the invention has been described above with reference to examples illus- 
trated in the attached drawings, it is to be understood that the invention is not limited 
thereto but can be modified in many ways within the scope of the inventive idea 
claimed in the following claims. 
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Claims 

1. A method for solving a system of linear equations, comprising the steps of: 

a. representing elements of a solution vector as fixed-point binary words 
each consisting of at least one bit; 

b. initialising the solution vector and an auxiliary vector; 

c. performing, for each bit representing the binary words, bit-wise itera- 
tions comprising the steps of: 

i. performing passes through all elements of the solution vector; 

ii. updating elements of the solution vector in the passes; 

iii. updating elements of the auxiliary vector in the passes; 

iv. repeating the passes until a finishing condition is fulfilled; 

d. stopping solving the system of linear equations when a stopping condi- 
tion is fulfilled. 

2. The method as defined in claim 1 wherein elements of the solution vector 
are initialised as zeros. 

3. The method as defined in claim 1 wherein the auxiliary vector is initialised 
by the right-side vector of the system of linear equations. 

4. The method as defined in claim 1 wherein the bit-wise iterations start from 
the most significant bit and proceed with the next less significant bit if the finishing 
condition is" fulfilled. 

5. The method as defined in claim 1 wherein in each pass, in turn, for each 
element of the solution vector a condition successful/unsuccessful is checked. 

6. The method as defined in claim 1 wherein the finishing condition is fulfilled 
if in a pass ho element of the solution vector is updated. 

7. The method as defined in claim 1 wherein the stopping condition is fulfilled 
if a predefined number of passes through all elements of the solution vector is exceeded. 

8. The method as defined in claim 1 wherein the stopping condition is fulfilled 
if a predefined number of bit-wise iterations, defining the number of valuable bits in 
elements of the solution vector as well as accuracy of the solution, is exceeded. 

9. The method as defined in claim 1 wherein the stopping condition is fulfilled 
if a computer time predefined for performing this method is finished. 
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10. The method as defined in claim 1 wherein when passing through elements 
of the solution vector the order of analysing elements of the solution vector in the pass 
is arbitrary. 

11. The method as defined in claim 1 wherein when passing through elements 
of the solution vector the pass starts from an element whose position corresponds to the 
position of an element of the auxiliary vector with maximum amplitude and in the order 
of reducing the amplitude. 

12. The method as defined in claim 11 wherein ordering elements of the auxil- 
iary vector is performed to define the order of elements in the pass. 

13. The method as defined in claim 5 wherein updating elements of the solution 
vector and the auxiliary vector is performed only if the condition success- 
ful/unsuccessful is successful. 

14. The method as defined in claim 13 wherein the only element of the solution 
vector is updated, for which the condition successful/unsuccessful is checked. 

15. The method as defined in claim 14 wherein a finite number of possible up- 
dates of the element of the solution vector are analysed for finding a preferable update 
and the element of the solution vector when updated is set to be equal to the preferable 
update. 

16. The method as defined in claim 15 wherein finding the preferable update 
comprises the steps of: 

e. calculating, for each possible update, an auxiliary value; 

f. finding a minimum among the auxiliary values; 

g. calculating a threshold; 

h. comparing the minimum with the threshold; 

i. choosing the preferable update as that corresponding to the minimum. 

17. The method as defined in claim 16 wherein the condition success- 
ful/unsuccessful is successful if the minimum among the auxiliary values is less than the 
threshold, and the condition successful/unsuccessful is unsuccessful if the minimum 
among the auxiliary values is higher than or equal to the threshold. 

18. The method as defined in claim 16 wherein calculating the auxiliary values 
is based on the corresponding element of the auxiliary vector. 
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19. The method as defined in claim 16 wherein calculating the threshold is per- 
formed by using a diagonal element of the coefficient matrix, the diagonal element cor- 
responding to the element of the solution vector, and a step-size parameter. 

20. The method as defined in claim 19 wherein the step-size parameter is de- 
creased after each bit- wise iteration. 

21. The method as defined in claim 20 wherein decreasing the step-size is by a 
factor of two. 

22. The method as defined in claim 13 wherein elements of the auxiliary vector 
are updated by using elements of the coefficient matrix and the step-size parameter. 

23. The method as defined in claim 22 wherein an element of the auxiliary vec- 
tor is updated by using elements of a row of the coefficient matrix, the row correspond- 
ing to the updated element of the auxiliary vector, and the step-size parameter. 

24. A computer system for solving a system of linear equations, comprising: 

j. a host processor producing itself or receiving from other devices pa- 
rameter signals representing elements of a coefficient matrix and right- 
side vector of the system of linear equations and transmitting to other 
devices parameter signals representing elements of a solution vector; 

k. a host bus coupled to the host processor; 

1. an internal bus; 

m. a first means for storing and updating elements of the solution vector, 
the first means coupled to the host bus and the internal bus; 

n. a second means for storing elements of the coefficient matrix, the sec- 
ond means coupled to the host and internal buses; 

o. a third means for determining successful iterations and preferable up- 
dates, the third means coupled to the internal bus and the second means; 

p. a fourth means for storing and updating elements of an auxiliary vector, 
the fourth means coupled to the host bus, the internal bus, and the sec- 
ond means. 

25. The computer system of claim 24 comprising a controller coupled to the 
host bus and the internal bus, receiving control signals from the host processor through 
the host bus and producing control signals for the internal bus. 
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26. The computer system of claim 24, wherein the first means contains a mem- 
ory means for storing elements of the solution vector and a means for adding or sub- 
tracting, and the first means updates elements of the solution vector by adding or sub- 
tracting a step-size parameter. 
5 27. The computer system of claim 24 5 wherein the second means contains a 

memory means for storing elements of the coefficient matrix and a means for bit- 
shifting, and the second means performs bit-shifts of elements of the coefficient matrix. 

28. The computer system of claims 24, wherein the fourth means contains a 
memory means for. storing elements of the auxiliary vector and a means for adding or 

10 subtracting, and the fourth means updates elements of the auxiliary vector by adding or 
subtracting bit-shifted elements of the coefficient matrix. 

29. The computer system of claim 24, wherein the fourth means receives ini- 
tialisation data from the host bus. 

30. The computer system of claim 24, wherein the fourth means receives ini- 
15 tialisation data from the host bus, the initialisation data being elements of the right-side 

vector of the system of linear equations. 

31. The computer system of claim 24, wherein the second means receives ele- 
ments of the coefficient matrix from the host bus. 

32. The computer system of claim 24, wherein the third means determines suc- 
- 20 cessful iterations by comparing elements of the auxiliary vector and bit-shifted elements 

of the coefficient matrix. 

33. A multiuser receiving method in a data transmission system in which code 
division multiple access, involving multiuser interference among respective signals, 
each signal representing a succession of data signals translated into bits and transmitted 

25 at a rate of a plurality of chips per bit, spread by a respective spreading code, is applied 
for detecting a particular data signal, from among a plurality of data signals, said 
method comprising: 

q. filtering matched with the spreading codes and applied to the received 
signal to obtain respective output signals; 
30 r. transforming the matched-filter output signals by solving a system of 

linear equations of the kind Rh=p where R is a NxN cross-correlation 
matrix of the spreading codes, p is a Nxl vector grouping the matched- 
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filter output signals, h is the Nxl solution vector representing the trans- 
formed signals, and N is a number of used spreading codes, wherein 
solving the system of linear equations comprises the steps of: 

i. representing elements of a solution vector as fixed-point binary 
words each comprising at least one bit; 

ii. initialising the solution vector and an auxiliary vector; 

iii. performing, for each bit representing the binary words, bit-wise 
iterations comprising the steps of: performing passes through all 
elements of the solution vector; updating elements of the solu- 
tion vector in the passes; updating elements of the auxiliary vec- 
tor in the passes; repeating the passes until a finishing condition 
is fulfilled; 

iv. stopping solving the system of linear equations when a stopping 
condition is fulfilled. 

s. subjecting the transformed signals to obtain estimates of the data sig- 
nals. 

34. A multiuser receiver in a data transmission system in which code division 
multiple access, involving multiuser interference among respective signals, each signal 
representing a succession of data signals translated into bits and transmitted at a rate of 
" a plurality of chips per bit, spread by a respective spreading code, is applied for detect- 
ing a particular data signal, from among a plurality of data signals, said multiuser re- 
ceiver comprising: 

t. filters matched with spreading codes contained in the received signals; 
u. a computer system for solving systems of linear equations of the kind 
Rh=p where R is a NxN cross-correlation matrix of the spreading codes, 
P is a Nxl vector grouping the matched-filter output signals, h is the 
Nxl solution vector representing the output signals of the computer sys- 
tem, and TV is a number of used spreading codes, wherein the computer 
system performs a sequence of operations comprising the steps of: 

i. representing elements of a solution vector as fixed-point binary 

words each consisting at least one bit; 
ii. initialising the solution vector and an auxiliary vector; 
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iii. performing, for each bit representing the binary words, bit-wise 
iterations comprising the steps of: performing passes through all 
elements of the solution vector; updating elements of the solu- 
tion vector in the passes; updating elements of the auxiliary vec- 
tor in the passes; repeating the passes until a finishing condition 
is fulfilled; 

iv. stopping solving the system of linear equations when a stopping 
condition is fulfilled. 

v. a means for estimating the data signals from the output signals of the 
computer system. 

35. A multiuser receiver according to claim 34, wherein the computer system 
for solving the system of linear equations comprises: 

w. a host processor receiving parameter signals representing elements of 
the cross-correlation matrix of the spreading codes and right-side vector 
of the system of linear equations and transmitting parameter signals rep- 
resenting elements of the solution vector; 

x. a host bus coupled to the host processor; 

y. an internal bus; 

z. a first means for storing and updating elements of the solution vector, 
the first means coupled to the host bus arid the internal bus; 

aa. a second means for storing elements of the cross-correlation matrix, the 
second means coupled to the host bus and the internal bus; 

bb. a third means for determining successful iterations and preferable up- 
dates, the third means coupled to the internal bus and the second means; 

cc. a fourth means for updating and storing elements of an auxiliary vector, 
the fourth means coupled to the host bus, the internal bus, and the sec- 
ond means. 

36. An adaptive filter for receiving a first signal from a first transmission line 
and a second signal from a second transmission line, said first signal partially leaking 
from said first transmission lineto said second transmission line as an echo, said adap- 
tive filter comprising a filter means coupled to said first transmission line and respon- 
sive to said first signal for producing an estimated echo signal determined in accordance 
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with filter coefficients, a coefficient generating means for generating said filter coeffi- 
cients, and subtracting means coupled to said filter means and connected in said second 
transmission line for subtracting said estimated echo signal from said second signal on 
said second transmission line so as to cancel said echo signal, said coefficient generat- 
ing means comprises: 

dd. a first means coupled to said first transmission line and responsive to 
said first signal for producing a series of autocorrelation coefficients of 
said first signal; 

. ee. a second means coupled to said first and said second transmission lines 
and responsive to said first and said second signal for producing a series 
of cross-correlation coefficients between said first signal and said sec- 
ond signal; and 

ff. a means coupled to said first and said second means for generating said 
filter coefficients from said autocorrelation and cross-correlation coeffi- 
cients to deliver said filter coefficients to said filter means, said third 
means is a computer system for solving a system of linear equations 
whose coefficient matrix comprises said autocorrelation coefficients and 
whose right-side vector comprises said cross-correlation coefficients, 
wherein said computer system for solving the system of linear equations 
20 performs a sequence of operations comprising the steps of: 

i. representing elements of a solution vector as fixed-point binary 
words each consisting of at least one bit; 

ii. initialising the solution vector and an auxiliary vector; 

iii. performing, for each bit representing the binary words, bit- wise 
iterations comprising the steps of: performing passes through all 
elements of the solution vector; updating elements of the solu- 
tion vector in the passes; updating elements of the auxiliary vec- 
tor in the passes; repeating the passes until a finishing condition 
is fulfilled; 

iv. stopping solving the system of linear equations when a stopping 
condition is fulfilled. 
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37. The adaptive filter of claim 36, wherein the computer system for solving the 
system of linear equations comprises: 

gg. a host processor receiving parameter signals representing elements of 
the coefficient matrix and the right-side vector and transmitting parame- 
ter signals representing elements of the solution vector; 

hh. a host bus coupled to the host processor; 

ii. an internal bus; 

jj. a first means for storing and updating elements of the solution vector, 
the first means coupled to the host bus and the internal bus; 

kk. a second means for storing elements of the coefficient matrix, the sec- 
ond means coupled to the host bus and the internal bus; 

11. a third means for determining successful iterations and preferable up- 
dates, the third means coupled to the internal bus and the second means; 

mm. a fourth means for storing and updating elements of an auxiliary 
vector, the fourth means coupled to the host bus, internal bus, and the 
second means. - 

38. The adaptive filter of claim 36, wherein said filter means is a transversal fil- 
ter. 
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ABSTRACT 

DATA PROCESSING, PARTICULARLY IN COMMUNICATIONS SYSTEMS 

The invention provides a method for solving systems of linear equations. Elements 
of a solution vector are represented as fixed-point binary words. The method comprises 
5 the steps of initialising a solution vector and an auxiliary vector and bit-wise iterations 
performing passes through elements of the solution vector, updating elements of the so- 
lution and auxiliary vectors in the passes, and repeating the passes until a required accu- 
racy of the solution is achieved. The method solves a system of linear equations by 
starting the bit- wise iterations with the most significant bit of the words and proceeds to 
10 a next (less significant) bit when the solution achieved for the current bit is not changed 
in the last pass. The solution is obtained with a computer system programmed to per- 
form a sequence of operations. 
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